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Abstract 

For a robust interpretation of upcoming observations from Planck and LHC 
experiments it is imperative to understand how the inflationary dynamics of 
a non-minimally coupled Higgs scalar field may affect the degeneracy of the 
inflationary observables. 

We constrain the inflationary observables and the Higgs boson mass dur- 
ing observable inflation by fitting the Hubble function and subsequently the 
Higgs inflationary potential directly to WMAP5+BAO+SNI We obtain for 
Higgs mass mmggs = 143.73„g 4 3i 7 GeV at 95% CL for the central value of top 
quark mass. We show that the inflation driven by a non-minimally coupled 
scalar field to Einstein gravity leads to significant changes of the inflationary 
parameters when compared with the similar constraints from the standard 
inflation. 

Key words: Key words: Ination, Higgs boson, Standard Model, 
Non-minimal coupling 
PACS: 98.80.Cq, 14.80. Bn 

1. Introduction 

The primary goal of particle cosmology is to obtain a concordant descrip- 
tion of the early evolution of the Universe consistent with both unified field 
theory and astrophysical and cosmological measurements. Inflation is the 
most simple and robust theory able to explain the astrophysical and cosmo- 
logical observations, providing at the same time self-consistent primordial 
initial conditions [l[ and mechanisms for the quantum generation of scalar 
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(curvature) and tensor (gravitational waves) perturbations 0]. 
In the simplest class of inflationary models, inflation is driven by a single 
scalar field <fi (inflaton) with some potential V(4>) minimally coupled to Ein- 
stein gravity. The perturbations are predicted to be adiabatic, nearly scale- 
invariant and Gaussian distributed, resulting in an effectively flat Universe. 
The WMAP 5-year CMB measurements alone |3| or complemented with other 
cosmological datasets j3] support the standard inflationary predictions of a 
nearly flat Universe with adiabatic initial density perturbations. In partic- 
ular, the detected anti-correlations between temperature and E-mode polar- 
ization anisotropy on degree scales [5| provide strong evidence for correlation 
on length scales beyond the Hubble radius. 

Alternatively one can look on inflationary dynamics based on the struc- 
ture of the Standard Model (SM) of particle physics. Recently, Bezrukov and 
Shaposhnikov [6| reported the possibility that the SM with an additional non- 
minimally coupled term of the Higgs field and Ricci scalar can give rise to 
inflation without the need for additional degrees of freedom already present 
in SM. The resultant Higgs inflaton effective potential in the inflationary 
domain allows for a range of cosmo logically values for the Hi ggs boson mass 
even after the inclusion of the radiative corrections 0, @, lilliol. 11]. More 



precisely, the present cosmological constraints on the Higgs mass are based 
on mapping between the Renormalization Group (RG) flow and the spectral 
index of the curvature perturbations emerging from the analysis of WMAP 
5-year data complemented with astrophysical distance measurements jij. 
However, after the publication of WMAP results, many works devoted to 
the reconstruction of the inflationary dynamics during observable inflation 



13l . \\M show the compatibility between WMAP results and the predictions 



of the standard inflation driven by a single scalar field with some potential 
V(4>) which is minimally coupled with Einstein gravity. 

In order to have a robust interpretation of upcoming observations from 
Planck 15J and LHC 16( experiments it is imperative to understand how 
the inflationary dynamics of a non-minimally coupled Higgs scalar field may 
affect the degeneracy of the inflationary observables. 

In this paper we aim to constrain the inflationary observables and the Higgs 
boson mass during observable inflation by fitting the Hubble function H(ip), 
and subsequently the Higgs inflationary potential V(<p), directly to WMAP 
5- year data 0, 0] complemented with geometric probes from the Type la su- 
pernovae (SN) distance-redshift relation and the baryon acoustic oscillations 
(BAO) measurements. 
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2. Perturbations driven by non- minimally coupled scalar field 



We briefly review the computation of the power spectra of scalar and 
tensor density perturbation generated during inflation driven by a single 
scalar field non-minimally coupled to gravity via the Ricci scalar R 17J, [18 , 

3, ET 
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Jordan frame. We start with the generalized action in the Jordan frame 17 



S 



d A x\ 



-g 



~U&)R-~G(tp)(VpT-V(tp) 



(1) 



where U((p) is a general coefficient of the Ricci scalar, G((p) is the general 
coefficient of kinetic energy and V((p) is the general potential. The general- 
ized U((p)H gravity theory in (pQ) includes diverse cases of coupling. 
For generally coupled scalar field U = (7 + £,mpf(p 2 ) and G(ip) = 1 where 
mpi ~ 2.4 x 10 18 GeV is the reduced Planck mass and 7 and £ are constants. 
The non-minimally coupled scalar field is the case with 7 = 1 while the con- 
formal coupled scalar field is the case with 7 = 1 and £ = —1/6. Throughout 
the paper we take sign(^) = — 1, keeping the same sign convention for £ as in 
Ref. |l7|. One should note that other authors 2l|, |22J use the opposite sign 
for £. The background equations in a flat Friedmann-Lemaitre-Robertson- 
Walker (FLRW) cosmology are given by: 



H 



HU-U 



Cp + 3Hp + — (G^ 2 - mj{(p\ v + 2 V{!p), v ) = 
2Cr 



(2) 

(3) 
(4) 



where a is the cosmological scale factor (a = 1 today), H is the Hubble 
expansion rate, overdots denotes the time derivatives and jip = d/dtp. 
Einstein frame. The conformal transformation to the Einstein frame for the 
action (CQ) can be achieved by defining the Einstein metric as fisj ] : 

9n,v = ^g^u , (5) 

where Q = U(ip) is the conformal factor. The action in the Einstein frame 
takes the canonical form: 



S, 



d x\/—g 



iR-i(Vy3) 2 -f(v,) 



(6) 
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where the new scalar field (p and the new potential V are defined through: 



3m 2 pi 



V u(<p) J 



+ 



G{tp) 



U{<p) 



V{<p) 



V{<p) 



<p=tp((p) 



U 2 (ip) 



The Higgs field suppression factor S{ip) defined in the above equation was 
first introduced in Ref. 23 1 in the context of CMB anisotropy generation in 
non-minimally inflation model to account for the suppression of the Higgs 
propagator due to non-minimal coupling with gravity. 

Decomposing the conformal factor into the background part and the per- 
turbed part: 

SQ(x,t) 



n(x,t) = n(t) 1 + 



Q(t) 



we get: 

a = aVQ , di = v^Qdt , H 




= ¥ + 



5Q 



(8) 



(9) 



In the above relations a, i, H, (p are respectively the cosmological scale 
factor, the cosmic time, the Hubble expansion rate and the scalar field in the 
Einstein frame. 

2.1. Cosmological perturbations in the Jordan frame 

Scalar perturbations (S). The equation of motion for the Lagrangian (jTJ 
is given by 181 ]: 



k 2 

a A Q s n) +—11 = 0, with Q s 



if + 3U 2 /2U 
(H + U/2U) 2 



(10) 



Here k is the comoving and 1Z is the intrinsic curvature perturbation of the 
co-moving hypersurfaces during inflation: 



where Sip is the perturbation of the scalar field (p. Neglecting the contribution 
of the decaying modes, the scale dependence of the amplitudes of the scalar 
perturbations can be obtained by integrating the mode equation jiil ]: 



<+ k 2 - 



u k = 



;i2) 



(7) 
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In the above equation the prime denote the derivative with respect to the 
conformal time 77(a) — J a(t)~ l dt. For scalar perturbations the term z" j z 
can be written as 1211: 



where 



(aHf 



[l + 5 s )(2 + 5 s + e) + ^- 

aid 



(13) 



H 
IP 



Qs 



2HQ, 



(14) 



Tensor perturbations (T). Tensor perturbations sa tisfy the same form of equa- 
tion as Eq. lflO!) with replacement Qs — > Qt — U 17]. The scale dependence 
of the amplitude of tensor (T) perturbations can also be obtained by inte- 
grating Eq. (ll2l) taking the term z" j z of the form 



z T 



{aHf 



;i + 5 T )(2 + 5 T + e) + ^E_ 



(15) 



where 



<h 



Q 



T 



2HQ T 

2.2. The equivalence of the perturbations spectra 
Introducing the following quantities: 



(16) 



dH 

dt 
dip 
H 



dQs/dt 
2HQ S 



Si 



dQnf j dt 
2HQ T 



(17) 



"-r \ _ Qs r v ( 

-) --jj, Qt~ — 



Qt 
U 



and making use of the relations given in Eq. fTB]) . it can shown that both scalar 
and tensor perturbations in the Jordan frame exactly coincide with those in 
the Einstein frame 25|, |26| implying that power spectra of the cosmological 
perturbations are invariant under the conformal transformation: 



V s {k) = V s (k) 



V T (k) = V T (k) . 



(18) 
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The numerical evaluation of the perturbations power spectra involves 
solving equations (fT21) . (ITS]) and (TTBT) for each value of the wavenumber k, 
the evolution of \v,k\/ 'zs,t to a constant value defining the observable power 
spectra Vs,t- Then the scalar and tensor spectral indexes ris,T and the run- 
ning of scalar spectral index as at the Hubble radius crossing k* are defined 
as: 



dlnV s (k) 



dink 



7ZT 



d\nP T (k) 



k t =aH 



dink 



dn s {k) 



k*=aH 



dink 



(19) 



k*=aH 



3. Higgs boson as inflaton 

Higgs field couplings in the inflationar 



add non-minimal coupling to gravity 22 



domain. Higgs boson as inflaton 
I, B 03 ■ Taking Higgs doublet 
in the form (l/\/2)(0, v + (p), the part of the general action given in Eq.([T]) 
relevant for inflation is only the scalar part. The quantum corrections due 
to the effects of particles of the Standard Model interacting with the Higgs 
boson through quantum loops modify the action coefficients U, G and V that 



in the Jordan frame take the form 10, 12 



m 



T „2t 



U(t) = l+Z(t)G 2 (t)^ T = l + at)G(t)^e 



m 



Pi 



m 



pi 



G(t) 
V(t) 



ex}) \ I dt'j(t' 
^p-G\t)(^ 2 (t) 2 -v 2 ) 2 



^G 4 (t)m*e 4i 



(20) 
(21) 
(22) 



where v ~ 246.2 GeV is the vacuum expectation value for the Higgs, A is the 
quadratic coupling constant of the Higss boson with the mass ran = V2Xv 
and j(t) is the Higgs field anomalous dimension. In the above equations the 
scaling variable 



t 



t(ip) = ln{ip/m T ) , <p(t) = m T e l (23) 
normalizes the Higgs field and all the running couplings to the top quark 



mass scale tut = 171.3 ± 1.1 ± 1.2 GeV [27|. As the energy scale of inflation 
is many order of magnitude above the electroweak scale (<p(t) » v), in the 
right hand side of Eq. fT22"]) the potential was approximated by a quadratic 
potential V((p) = \(p 4 /4: . 
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The Higgs field suppresion factor S(t), the Higgs potential V(t) and the first 
two potential slow-roll (PSR) parameters ey and r\y in the Einstein frame 
are given by: 



G(t)U(t)+3U^ t) ' w XP{t) 



ey^m^i^j S-\ Vv = 2e + m Pl 6 j^S-\ (25) 

The amplitude of the scalar (curvature) perturbations at the Hubble radius 
crossing fc* is related to the Higgs potential in the Einstein frame through: 



A 



2 = V_ 

s 24vr 2 - 4 



nr pl ey 



(26) 

fc. 



We compute the various t-dependent running constants, the Higgs field prop- 
agator suppretion factor and Higgs field anomalous dimension by integrating 
the two-loop renormalization group /5-functions for SU{2) x £7(1) gauge cou- 
plings g',g, the SU(3) strong coupling g s , the top Yukawa cou plin g h t and 
the Higgs quadratic coupling A, as collected in the Appendix of [291 ] . For the 
non-minimal coupling constant £ we use the /3-function from the Appendix 
of M. At t = 0, which corresponds to the top quark mass scale, the top 



Yukawa coupling h t (0) and Higgs quadratic coupling A(0) are fixed by their 
pole masses [29]. The gauge coupling constants at itlt scale are taken from 
@| : # 2 (0) = 0.4202, g' 2 (0) = 0.1291 and # s 2 (o) = 0.1.3460. The value of £(0) 
is determined so that at the initial point of the slow roll inflation, ti, the 
non-minimal coupling constant £(U) is such that the calculated value of the 
amplitude of density perturbationsgiven in Eq. fl26l) agrees with the measured 
result obtained by WMAP team [4|. 

Higgs field equations of motion. Hereafter we will work in the Einstein 
frame dropping the 'hat' when express different quantities in this frame. The 
inflationary dynamics and the cosmological perturbations power spectra are 
entirely determined by the dynamics of the Higgs field <p(t). 
By using the transformation relations from Jordan frame to Einstein frame 
given in Eq.Q, to the leading order in slow- roll approximation (ip « 3H(f, 
neglecting also the higher powers in (p ) the general background Eqs. (|2D - 



7 



flU) reduce to |20j,|2l|:. 



a J 
UV 



2 



1 




H 2 



3ml 



(27) 



PI L 



3Hip 




-2 



(28) 



It follows that the amount of the expansion from any epoch to the end of the 
inflation (ey = 1) that is controlled by the number of e-folds J\f is: 



The CMB angular power spectra. We obtain the CMB temperature 
anisotropy and polarization power spectra by integrating the mode equa- 
tion (TT2"j) and the Higgs field equations (|2"T1) and (I2"g|) with respect to the 
conformal time. We consider wavenumbers in the range [5 x 10~ 6 — 5] Mpc -1 
needed to numerically derive the CMB angular power spectra and the Hub- 
ble crossing scale fc* = 0.002Mpc _1 . At this scale, the value of the Higgs 
scalar field (p* corresponds to Af = k*/aH ~ 59 e-folds till the end of the 
inflation and the value of the scaling parameter is ~37 0, S], as indicated 
by the central value of the amplitude of scalar perturbations measured by 
WMAP at k* [3]. For each wavenumber k in the above range our code inte- 
grates Eqs. f|T2|) . (T27j) and ff28|) and the (3 functions of the t-dependent running 
constant couplings in the observational inflationary window imposing that k 
grows monotonically to the wavenumber k*, eliminating at the same time the 
models violating the condition for inflation < e# = —H/H 2 < 1. 
We reconstruct the Hubble expansion rate H (0 — 0*) from the data by using 
the Taylor expansion up to the cubic term, equivalent to keeping the first 
three slow-roll parameters. 

For comparison, we compute the CMB angular power spectra in the standard 
inflation model with minimally coupled scalar field with a general potential 
V((p). The details of this computation can be found in 




(29) 
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Figure 1: ID marginalized likelihood probability distributions for seven independent pa- 
rameters obtained from the fit to WMAP5+SN+BAO dataset of the inflation model with 
non-minimally coupled Higgs scalar field (continuous lines) and of the standard infla- 
tion model with minimally coupled scalar field (dashed lines). The last plot shows the 
value of Hubble parameter. All parameters are computed at the Hubble radius crossing 
/c*=0.002Mpc- 1 . 



4. Results 

We use the Markov Chain Monte Carlo (MCMC) technique to recon- 
struct the Higgs field potential and to derive constraints on the inflationary 
observables and the Higgs mass from the following datasets. 
The WMAP 5-year data [4| complemented with geometric probes from the 
Type la supernovae (SN) distance-redshift relation and the baryon acoustic 
oscillations (BAO). The SN distance-redshift relation has been studied in 
detail in the recent unified analysis of the published heterogeneous SN data 
sets - the Union Compilation08 [32] . The BAO in the distribution of galaxies 
are extracted from the Sloan Digital Sky Surveys (SDSS) and Two Degree 
Field Galaxy Redshidt Survey (2DFGRS) Q. The CMB, SN and BAO data 
(WMAP5+SN+BAO) are combined by multiplying the likelihoods. 
We decided to use these measurements especially because we are testing 
models deviating from the standard Friedmann expansion. These datasets 
properly enables us to account for any shift of the CMB angular diameter 
distance and of the expansion rate of the Universe. 
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Figure 2: The figure presents the results of the fit to WMAP5+SN+BAO dataset of the 
inflation model with non-minimally coupled Higgs field in blue (continuous lines) and 
the results of the fit to the same dataset of the standard inflation model with minimally 
coupled scalar field in yellow (dashed lines). The top plot in each column shows the 
probability distribution of different scalar inflationary observables while the other plots 
show their joint 68% and 95% confidence intervals. All parameters are computed at the 
Hubble radius crossing fc*=0.002Mpc _1 . 
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Fi gure 3: 2D marginalized likelihood probability distributions in the planes TTiHiggs 

- £ (left 

panel) and mHi ggs -\ (right panel) obtained from the fit to WMAP5+SN+BAO dataset 
of the inflation model with non-minimally coupled Higgs scalar field. 



The likelihood probabilities are evaluated by using the public packages 



CosmoMC and CAMB |33|, |34| modified to include the formalism for in- 
flation driven by non-minimally coupled Higgs scalar field, described in the 
previous sections. Our fiducial model is the ACDM standard cosmological 
model described by the following set of parameters receiving uniform priors: 

H' 2 H" H"'H' 



n b h 2 



A 2 



H 2 ' jj ' R 2 

where: Qbh 2 is the physical baryon density, fl c dmh 2 is the physical dark matter 
density, 9 S is the ratio of the sound horizon distance to the angular diameter 
distance, r is the reionization optical depth, m Higgs is the Higgs boson mass 
and H', H" and H'" are the derivatives of the Hubble expansion rate with 
respect to the Higgs scalar field. 

For comparison, we use the MCMC technique to reconstruct the standard 
inflation field potential and to derive constraints on the inflationary observ- 
ables from the fit to WMAP5+SN+BAO dataset of the standard inflation 
model with minimally coupled scalar field and a general potential V(<f>) 
For this case we use the same set of input parameters with uniform priors as 
in the case of non-minimally coupled Higgs scalar field inflation, except for 
Higgs mass. 

For each inflation model we run 120 Monte Carlo chains, imposing for each 
case the Gelman & Rubin convergence criterion fiF 
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Fig. 1 presents the ID marginalized likelihood probability distributions of 
the fiducial model parameters obtained from the fit to WMAP5+SN+BAO 
dataset of the inflation model with non-minimally coupled Higgs scalar field 
compared with the similar distributions obtained from the fit to same dataset 
of the standard inflation model with minimally coupled scalar field. One can 
see the differences between the likelihood probabilities caused by the differ- 
ence in the dynamics of the scalar fields during inflation. One should note the 
smaller value of the Hubble expation rate during inflation for the inflation 
model involving the Higgs field as inflaton. 

In Fig. 2 we present the constraints on the inflationary parameters, namely 
the scalar spectral index ns, the running of the scalar spectral index as, the 
amplitude of the scalar (curvature) density perturbations A 2 S and the ratio of 
tensor to scalar amplitudes R, as obtained from the fit to WMAP5+SN+BAO 
dataset of the inflation model with non-minimally coupled Higgs field com- 
pared with the constrains obtained from the fit to the same dataset of the 
standard inflation model with minimally coupled scalar field. We also show 
the constraints obtained on the Higgs mass from our analysis. The mean val- 
ues and the 95% CL intervals of these parameters are presented in Table 1. 
We obtain a Higgs mass value of mmggs = 143.73_g 4 3^ 7 GeV at 95% CL for 
the central value of top quark mass. 

Our result is compatible with the previous constrains of the Higgs mass from 
cosmological data 0, S, @, 0, 11] however the upper and lower limits of the 



Higgs mass obtained from this analysis are tightly constrained. 
Table 1 clearly show the differences between the estimates of the inflationary 
parameters obtained from the fits of the inflation models. 
Fig 3 shows the 2D marginalized likelihood probability distributions in the 
planes muiggs - £ and m Higgs -X. One should note the strong correlation 
between Higgs mass and the Higgs quadratic coupling A and the larger de- 
generacy between the Higgs mass and non-minimal coupling constant £, in- 
dicating the existence of a degeneracy between £ and A that affect the Higgs 
mass value inferred from cosmological and astrophysical measurements. 

5. Conclusions 

In order to have a robust interpretation of upcoming observations from 



Planck [15| and LHC [16( experiments it is imperative to understand how 
the inflationary dynamics of a non-minimally coupled Higgs scalar field may 
affect the degeneracy of the inflationary observables. 
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We constrain the inflationary observables and the Higgs boson mass during 
observable inflation by fitting the Hubble function H(ip), and subsequently 
the Higgs inflationary potential V((p), directly to WMAP 5-year data com- 
plemented with geometric probes from the Type la supernovae (SN) distance- 
redshift relation and the baryon acoustic oscillations (BAO) measurements. 
We obtain a Higgs mass value of m Higgs = 143.73_g 4 3| 7 GeV at 95% CL for 
the central value of top quark mass mT op =171.3 GeV. Our result is compat- 
ible with the previous constrains of the Higgs mass from cosmological data 
@, H, @, 0, however the upper and lower bounds of the Higgs mass 



obtained from our analysis are tightly constrained. 

The strong correlation between Higgs mass and the Higgs quadratic cou- 
pling A and the larger degeneracy between the Higgs mass and non-minimal 
coupling constant £, indicate the existence of a degeneracy between £ and A 
that affect the Higgs mass value inferred from cosmological and astrophysical 
measurements. 

We also show that the inflation driven by a non-minimally coupled scalar 
field to Einstein gravity leads to significant changes of the inflationary param- 
eters when compared with the similar constraints from the standard inflation 
with minimally coupled scalar field. 
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Table 1: The mean values and 95% CL lower and upper intervals of the derived pa- 
rameters obtained from the fit of the inflation model with non-minimally coupled Higgs 
scalar field and the standard inflation model with minimally coupled scalar field to 
WMAP5+SN+BAO dataset. All parameters are computed at the Hubble radius crossing 
/c»=0.002 Mpc" 1 . 



Model 


Standard Inflation Higgs Inflation 


Parameter 





fl b h 2 


n nooO.023 

u.uzz 021 


nooO.024 
U.UZOg 022 


n c h 2 


O-lllo'ioe 7 




T 


0-082^1 


D DQ4 0135 
u.u»^ 034 


0s 


±.uoy 1 034 


1 n44 L04S 


ln[10 10 A|] 


3-143 3 ; 083 


o oi 7 3.252 
°- 0± ' 2.999 


n s 


0.95683™ 


n 0770.997 


a s 


_n mo 0.021 
U.U1Z_ 03g 


n nno 0.0133 
U.UUO_ .0021 


R 


< 0.556 


< 0.234 



m H iggs(GeV) 


1 AO 7Q7158.707 


A 


o.i688il§ 


£ x lO" 5 


i-87i8:ffg 
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